T				! Darcy (T for Darcy-2D, F for FEM-2D)
0	1	0	1	! Xmin, Xmax, Ymin, Ymax
20	20			! Nx, Ny, (Ns=1)
4.	0.	0.	1.	! Dxx, Dxy, Dyx, Dyy
2	2			! Nxx, Nyy for f(x), g(y)
0	1			! coordinate x for f(x)	(xx(k), k=1,Nxx)
0	0			! f(x) value, s=1=Ns (fx(s,k), k=1,Nxx)
0	1			! coordinate y for g(y)	(yy(k), k=1,Nyy)
0	0			! g(y) value, s=1=Ns (gy(s,k), k=1,Nyy)
F	0	0		! delta (T or false) and its Xdelta, Ydelta
0				! C(1) IC (cons for all nodes)
'D'	'N'	'D'	'D'	! (BCt(f), f=1,4) (south, east, north, west)
21	2	2	2	! (Nxy(f), f=1,4) (south, east, north, west)
0	0.05	0.1	0.15	0.2	0.25	0.3	0.35	0.4	0.45
0.5	0.55	0.6	0.65	0.7	0.75	0.8	0.85	0.9	0.95
1				! (xy(1,i)  , i=1,Nxy(1))
0	0.078459096	0.156434465	0.233445364	0.309016994	0.382683432
0.4539905	0.522498565	0.587785252	0.649448048	0.707106781
0.760405966	0.809016994	0.852640164	0.891006524	0.923879533
0.951056516	0.97236992	0.987688341	0.996917334
1				! (vl(1,s,i), i=1,Nxy(1), s=1=Ns)
0	1			! (xy(2,i)  , i=1,Nxy(2))	
0	0			! (vl(2,s,i), i=1,Nxy(2), s=1=Ns)	
0	1			! (xy(3,i)  , i=1,Nxy(3))	
0	0			! (vl(3,s,i), i=1,Nxy(3), s=1=Ns)	
0	1			! (xy(4,i)  , i=1,Nxy(4))	
0	0			! (vl(4,s,i), i=1,Nxy(4), s=1=Ns)	
--------------------------------------------------------------------------------
! Conduction in an orthotropic rectangle												
! from	H.S. Carlslaw and J.C. Jaeger
! Conduction of Heat in Solids, 2nd. ed., Oxford, 1959, pp. 166-169											
! case V	(Eq. 26) p. 169											
! with specific choice: f(x) = sin (m pi x/a) so that Am=1 for m=n and Am=0 for m#n
! 0 < x < a/2=1, 0 < y < b=1, Ne= 20x20 (dx=dy=0.05)
! Dxx = k^2*Dyy, k=2
! generated by MeshGenRect v54.3 of 22/08/18
--------------------------------------------------------------------------------
formers:
--------------------------------------------------------------------------------

